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Late time Cosmological Phase transition and Galactic Halo as Bose-Liquid 



(N ■ 
Os ■ 

o\ : 



oo 



> 

00 

o 

(N 

in 
o 

On 



i 

Oh! 
CD . 
43 ■ 



X 



Sang- Jin Sin' 

Department of Physics, University of Florida, Gainesville FL. 32611 

(April 1992) 

We consider the ultra light pseudo Nambu-Goldstone boson appearing in the late time cosmolog- 
ical phase transition theories as a dark matter candidate. Since it is almost massless, its nature is 
more wave like than particle like. Hence we apply quantum mechanics to study how they form the 
galactic halos. Three predictions are made; (l)the mass profile p ~ r -1 ' 6 , (2)there are ripple-like 
fine structures in rotation curve, (3) the rotation velocity times ripple's wave length is largely galaxy 
independent. We compare the rotation curves predicted by our theory with the data observed. 



1. Introduction - Recently, motivated by the large 
scale structure, an ultra light Nambu-Goldstone boson 
was introduced in late time cosmological phase transi- 
tion theory. [1] The bottom line of the theory is that if a 
phase transition happens after decoupling, one can avoid 
the constraint imposed by isotropy of microwave back 
ground. If an event happens so late then the universe 
is already big, hence the relevant particle whose comp- 
ton wave length provides the scale of interesting struc- 
ture must be very light and this particle could be a dark 
matter candidate. [2] 

The dark matter distribution is most clearly imprinted 
in the rotation curve. The overall 'flatness' of the rota- 
tion curve (RC) is equivalent to mass profile p ~ r~ 2 , 
which can be obtained by assuming a thermal distribu- 
tion. [3] [4] Another way suggested to understand the 
dark matter distribution is a spherical infall model [5] [6] 
[7] [8] [9] , where the mass profile is given by p ~ r - 2 - 25 j 
giving a slowly decreasing rotation curve. However, the 
observed rotation curves are actually slightly increasing 
[10] though they are said to be 'flat'. So neither the 
thermal distribution nor the infall model is completely 
satisfactory. 

If dark matter consists of these particles whose mass 
is, say, m ~ 10 _24 eV r and dark matter density around 
us is ~ 10 -25 g/cm 3 , then the inter-particle distance is 
order of 10 -13 cto, while the compton wave length is of 
lOpc order. Hence there is no point talking about the 
individual particle's position and momentum. The dark 
matter is more like a wave than a particle, and the galac- 
tic halos are giant systems of condensed bose liquid. Here 
we are in the region reigned by quantum mechanics. In 
this context, we will consider the dark matter distribu- 
tion quantum mechanically and this is the purpose of our 
paper. Our result will show mass profile p ~ r -1 - 6 lead- 
ing to slightly increasing rotation curve. Further more 
due to the system's wave nature, rotation curve has the 
ripple-like fine structure. By the uncertainty principle, 
our theory predicts that the rotation velocity times the 
ripple's wave length is galaxy independent. 



sider the spherically symmetric case for simplicity and 
also assume that there are no inter-particle interactions 
apart from gravity. The classical energy of a particle 
moving inside a potential generated by other particles is 
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where M(r) = dr'Anr' 2 p(r'). Quantum mechanically, 
we are interested in a wave equation which governs the 
collective behavior of highly correlated bosonic particles. 
In principle, one has to work out Hartree's self consistent 
procedure, but here the fact that the system is a highly 
correlated system of bosons simplifies the life. Assuming 
that all or significant fraction of particles are in the same 
state, we take the condensation wave function ip [11] as 
the Hartree's self-consistent field. Then the density of 
particles is proportional to \ip\ 2 and we obtain a wave 
equation given by 
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2. Quantum considerations of dark matter 



We con- 



Aip+GmM — dr" ^r" 2 \i>\ 2 -i>{r) 
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(2) 

Here Mo in M — Mq J dr4:Trr 2 \ip\ 2 is a mass parameter 
introduced for dimensional reasons. The equation is non- 
linear, hence we can not superpose two solutions, nor do 
we have freedom to normalize it. However, as we shall 
show later, this equation has a scaling symmetry which 
allows us to resolve the ambiguity due to the choice of 
Mq, That is, the physical quantity M does not depend 
on the choice of Mq, which is a nontrivial property. 

Another remark is on the non-relativistic treatment. 
(From now on we set % = c = 1 except for emphasis.) 
One might wonder that when the mass of the particle 
is so small, whether a non-relativistic treatment is cor- 
rect. The justification is two-fold. The basic reasons are 
(l)the virial theorem tells us v 2 = GM/r is independent 
of individual particle mass m and (2)the total mass of the 
system is small enough as we shall show below. When 
the effective coupling constant e 2 := GmM ~ 1 or equiv- 
alently when M ~ l/Gm, the gravitational Bohr radius 
ro ~ 1/me 2 is equal to the compton wave length 1/m, 
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and the size of the system is equal to the size of a parti- 
cle. This means that particles lie on top of one another, 
which is nothing but a black hole system. In fact, the 
Schwarzschild radius is 2GM <~ 1/m, confirming our pic- 
ture. When this happens, the non-relativistic treatment 
breaks down. In our case, however, 1/Gm <~ 10 14 M Q 
while M ~ 1O 12 M so that v ~ e 2 ~ 1CT 2 . Hence the 
non-relativistic treatment is well justified. 

We are interested in a stationary solution ip(r,t) = 
e~ lEt / h tp{r). After scaling 



r = r r, ijj = r 3/2 tp, E = 



2m 1 



(3) 



where r = h 2 /2GM m 2 , we can rewrite (2) in terms of 
the radial wave function u(r), 



i"(f) + (e-/ df'^l df"u 2 (f"))u(f)=Q, 
Jo r Jo 



(4) 



where 4>{f) = -^=^p-. We have solved this equation nu- 
merically. Figures l(a)-(b) are the stationary solutions 
and corresponding velocity curves v(r) — v(r)/v — 
(j Q f u 2 df' /f) 1 / 2 where v a = ^/GM /r . 

Figure 1 (a),(b) 

We are interested in a higher state with nodes rather 
than the ground state, because there is no way for the 
system to lose its energy to relax to its ground state, 
since there is no dissipation mechanism. Let's define 
the node number n as the number of nodes plus 1. For n 
greater than 4, the rotation curve resembles the observed 
ones and the characteristic features do not change as we 
increase the number of nodes. The numerical study show 
that mass profile p <~ y-- 1 - 56 ^ >> \ f or n — g, and the 
power depends on n very weakly. 

We can 'understand' (not derive) behaviors qualita- 
tively by simple semi-classical analysis [12] Since M(r) ~ 
rv 2 , flatness of rotation curve is equivalent to M'(r) <~ u 2 
being constant. Usual WKB method gives 
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where p r = ^/2m(e — V(u)). Since the potential is in- 
creasing very slowly as a function of r, so is 1/^/p^; hence 
it gives a minor modification to the 'over-all flat' cosine 
curve. This explains the overall flat but slowly increasing 
aspect of the rotation curves. Also notice that the last 
wave is almost as long as the sum of all the others since 
p r near the 'turning point' is singular. The ripple like 
structure is an inevitable phenomena predicted in this 
approach. It is a direct consequence of using the wave 
theory under the one state (or small number of states) 
dominance asumption. 



Apart from mass profile and ripple's existence, our 
theory has one quantitative prediction which could prove 
or disprove it with more precise measurements of rota- 
tion curves. 

Rotation velocity times the ripple wave length is largely 
Galaxy independent. 

This is nothing but re-phrasing the uncertainty principle. 
More precisely, 



Vrotation^ripple — r— f\f)i 
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where /(f) is a slowly increasing dimension-less function 
given by 



m = 
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where V(f) = V(r)/mv 2 and r t is r t is the effective 
classical turning point. If we further assume that the 
potential changes very slowly as we change n » 1, we 
get the eigenvalue distribution e n <~ eo + ( n + /3) 2Q /( 2+a ) 
for V ~ r a . Numerical study shows that without includ- 
ing visible matter, a — 0.44 for n=6, and 0.45 for n=5. 
When we compare with data we have to include visible 
matter and this leads toa~ 0.39. See figure 3(c). 

Figure 2 (a),(b) 

Next we observe that as we vary the initial condition 
t/jo, the scale free quantities e and M := J u 2 dr have a 
simple dependence on it. Namely 

M = 20.3068^^, e = 1.4900^ for n = 5 (8) 

The coefficients depend on n, e.g, for n= 6 (5 nodes), 
they are given by 24.4408 and 1.5385 respectively. To 
understand this we have proved that equation (2) has 
the following symmetry. 

Scaling symmetry : if ^>(r, e) is a solution, so is 
A 2 V(Ar,A 2 e) 

The relation (5) is just a corollary of this rather interest- 
ing property. Also, by choosing A such that the norm of 
ip is equal to one, one can equate Mq and M. This is not 
convenient in practical numerical integration, however. 
We will discuss this property when we discuss the data 
fitting. 

3. Comparison with observations — We now compare 
the the rotation curves predicted in our theory and the 
actual observations. First we have to include the visible 
matter's effect. See figure 3(a), (b). Equation (2) has to 
be modified to be 



-1 f r dr' 

id t 4> = - — Aib + Gm 
2m 



dr"4n(Mo\i>\ 2 +p vis )-^ 
(9) 
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We choose the effective density of visible matter as the 
Plumber's potential for the bulge and a Yukawa type 
exponential decay for the disk. By including the visible 
matter's contribution, the rotation curve becomes more 
flat, but the wavelength of the ripple docs not change. 



Figure 3 (a),(b),(c),(d) 



The ripple structure is manifest for relatively few 
galaxies. The data of Rubin et.al showed that in any 
galaxy there are fluctuations that could indicate the fine 
structure, but most of the data is plagued with error 
bars that wash out the signal of the fine structure. The 
data is good enough for the rotation velocity's magnitude 
but too crude for the fine structure. Among 63 galaxies 
listed in [10], there were about five galaxies which have 
relatively small error bars and show ripples. These are 
NGC2998, NGC1357, NGC4594, NGC1620, NGC801. 
Galaxy NGC2998's data is particularly interesting for 
us and we plot the data and the theoretical curve in 
Figure 3(d). There are four ripples in the NGC2998's 
data but we suspect that out side the measured region 
there can still be density peaks and nodes that could 
lead to further ripples. However as n increases, the total 
galaxy mass increases; and for large enough n, f2 exceeds 
1, violating bias of present theoretical community. For 
example, if n is larger than 7, the total mass of the halo 
is bigger than 1O 13 M0. Hence there is some restriction 
on n. Here we take the minimal choice n=5 (four nodes). 

By measuring the ripple wavelength and the rotation 
velocity, we can deduce the total mass M of the Galaxy 
and the constituent particle mass m. For convenience we 
choose ipo — 0.2. Comparing the velocity and wavelength 
of, say, the second ripple, we have two relations, 

7r = 8kpc, and 0.32w = 204km/sec. (10) 

Solving these two, we get 

to = 3.3 x W- 23 eV, and M = 0.69 x 10 12 M Q . (11) 

Using the numerically calculated value M = 8.73575, 

M = M ■ M = 5.9 x 1O 12 M (12) 

The central density of dark matter is given by 

Pc = l/4^ | 2 (M /r 3 ) = 1.08 x l$- 22 g/cm\ (13) 

Therefore we can weigh the Galaxy by looking at the 
rotation curve only! In the model adopted, the visible 
matter's total mass is 2.7 x IQ^Mq. For other galaxies, 
by similar analysis the prediction of eq. (8) is roughly 
satisfied. More massive galaxies have larger rotational 
velocity hence shorter ripple wavelength. For future data 
analysis we give the general formula for m, M with n=5. 

8kpc204km/sec _ , t , 
to = — f x 2.7 x 10~ 23 eU (14) 



M = — - ( — — - . — f x 5.9 x 1O 12 M . (15) 
8kpc 204km/ sec' U V ; 

Now we show that above quantities are independent 
of the choice of the initial condition ipo or equivalcntly 
independent of the choice of M . Recall that under the 
scaling, ip — > X 2 tp, f — > f/A and e — > A 2 e must follow. 
Therefore had we chosen A^o instead of i/> above, we 
would get 7kpc/\ — 8km/ sec and 0.32w A = 200fcTO/sec, 
so that we get the same to and 1/A times M above. 
On the other hand M = J d 3 f\tp\ 2 must scale by A 4 /A 3 , 
therefore M — M ■ Mo is invariant under the scaling. 
The same argument shows p c is scale free. Hence we can 
choose large to save computing time. Lastly, had we 
chosen n larger than 7, total mass for halo is larger than 
10 13 M Q . 

4- Conclusion — In this paper we investigated the con- 
sequence of abundant existence of Nambu-goldstone bo- 
son appearing in late time phase transition. Since the 
particle wavelength is far larger than the inter-particle 
distances, galactic halos made of this species are highly 
correlated bose liquid. We set up a version of non-linear 
shrodinger equation describing the collective behavior of 
the system. Three predictions are made; (l)mass pro- 
file p ~ r~ 16 , (2)there are ripple-like fine structures in 
rotation curve, (3) rotation velocity times ripple's wave 
length is largely galaxy independent. We compare the 
rotation curves predicted by our theory with the data 
observed. We fine-tuned the mass of the particle to fit 
the data for halo and determine the total mass of the 
galaxy NGC2998. Among 62 rotation curve data of Ru- 
bin et.al., there is a one curve where all the details of 
the predictions are realized, but most of other data is 
not accurate enough to definitely confirm the the theory. 
We suggests that more precise measurement of rotation 
curve is highly worth-while to be done. If future obser- 
vation shows the ripple structures in many galaxies in 
the way predicted in this paper, it would be a strong 
indication that late-time phase transition is real. 
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